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The problem of how many entangled or, respectively, sep- 
arable states there are in the set of all quantum states is in- 
vestigated. We study to what extent the choice of a measure 
in the space of density matrices g describing iV-dimensional 
quantum systems affects the results obtained. We demon- 
strate that the link between the purity of the mixed states 
and the probability of entanglement is not sensitive to the 
measure chosen. Since the criterion of partial transposition 
is not sufficient to distinguish all separable states for N > 8, 
we develop an efficient algorithm to calculate numerically the 
entanglement of formation of a given mixed quantum state, 
which allows us to compute the volume of separable states for 
N — 8 and to estimate the volume of the bound entangled 
states in this case. 

03.65.Bz, 42.50.Dv, 89.70,+c 



I. INTRODUCTION 

Entangled states have been known almost from the 
very beginning of quantum mechanics and their some- 
what unusual features have been investigated for many 
years. However, the recent developments of the theory 
of quantum information and quantum computing caused 
a rapid increase of the interest in studying their prop- 
erties and possible applications. To illustrate this trend 
let us quote some data from the Los Alamos quantum 
physics archives. In 1994 only one paper posted in these 
archives contained the key word 'entangled' (or entan- 
glement) in the title, while two such papers were posted 
in 1995. Since then the number of such papers has in- 
creased dramatically, and was equal to 8, 30 and 70, in 
the consecutive years 1996, 1997 and 1998, respectively. 

We do not dare to fit some fast growing curves to these 
data nor to speculate, when such an increase will eventu- 
ally saturate. On the other hand, since so many authors 
have dealt with entangled states, it is legitimate to ask, 
whether such states are 'typical' in quantum theory, or if 
they are rather rare and unusual. Vaguely speaking, we 
shall be interested in the relative likelihood of encoun- 
tering an entangled state jy. One may also ask a com- 
plementary question concerning the set of the separable 
states, which can be represented as a sum of product 
states. 

Consider a quantum system described by the density 
matrix p which represents a mixture of the pure states of 
the N dimensional Hilbert space. Let us assume that the 



system consists of two subsystems, of dimension tia and 
Ub, where N — tiatib. To formulate the basic question 
'What is the probability of finding an entangled state of 
size N?', one needs to: 

(i) define the probability measure fi, according to 
which the random density matrices p are drawn, 

(ii) find an efficient technique, which would allow one 
to judge, whether a given mixed state p is entangled. 

Representing any density matrix in the diagonal form 
p = UdU^ we proposed M to use a product measure 
p = Ai x v, where Ai describes the uniform measure on 
the simplex Xa=i ^» = 1 and v stands for the Haar mea- 
sure in the space of unitary matrices U(N). Based on the 
partial transposition criterion [g,[|] we found that under 
this measure the volume of the set of separable states is 
positive and decreases with the system size N jjj . Some 
more general analytical bounds were also provided by Vi- 
dal and Tarrach 0. Recently Slater suggested estimat- 
ing the same quantity using some other measures in the 
space of the density matrices [plpl . One may thus expect 
that the volume of the separable states depends on the 
measure chosen. We show that this is indeed the case. In 
this work we investigate which statistical properties de- 
scribing the set of the entangled states may be universal, 
e.g. do not depend on the measure used. In particu- 
lar we demonstrate that the relation between the purity 
of mixed states and the probability of entanglement is 
not very sensitive to the measure assumed. Basing on 
numerical results we conjecture that the volume of the 
separable states decreases exponentially with the system 
size N. 

For N = 4 and JV=6a density matrix is separable if 
and only if its partial transpose is positive 13]. For N > 6, 
however, there exist states, which are not separable and 
do satisfy this criterion [f7J. These states cannot be dis- 
tilled to the singlet form and are called bound entangled 
P-|ll| . Since there are no explicit conditions allowing one 
to distinguish between separable and bound entangled 
states, in Rcf. jl| only the upper bound for the volume of 
separable states was considered for N > 6. In this paper 
we present an efficient numerical method of computing 
the entanglement of formation E [ [12] for any density ma- 
trix. This method allows us to estimate the volume of 
bound entangled states, by taking a reasonably small cut- 
off entanglement E c and counting these states satisfying 
the partial transposition criterion for which E > E c . Our 
numerical results are to large extent independent on the 
exact value of E c . 



The paper is organized as follows. In section II we 
review the necessary definitions and study how the upper 
bound of the volume of the separable states depends on 
the system size and the measure used. The next section is 
devoted to analysis of the simplest case N = 4, for which 
the bound entangled states do not exist. In this case the 
analytical formula for the entanglement of formation is 
known |13|Jl5| and we study how this quantity changes 
with the purity of the mixed states. In section IV we 
study the case N = 8 and estimate the volume of the free 
entangled states, bound entangled states and separable 
states. The paper is concluded by Section V containing 
a list of open questions. In Appendix A we prove the 
rotational invariance of the two distinguished measures 
A and A u , defined on the TV — 1 dimensional simplex 
and demonstrate the link with the ensembles of random 
matrices. The algorithm of computing the entanglement 
of formation for a given density matrix is presented in 
Appendix B. 



II. VOLUME OF THE STATES WITH POSITIVE 
PARTIAL TRANSPOSITION 

A. Product measures in the space of mixed density 
matrices 

To discuss the probability, that a mixed quantum state 
possesses a given property, one needs to define a proba- 
bility measure \x in the space of density matrices (N x N 
positive Hermitian matrices with trace equal to unity). 
Each density matrix can be diagonalized by a unitary 
rotation. Let B be a diagonal unitary matrix. Since 
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UdU^ = UBdB ] U\ 



(1) 



therefore the rotation matrix U is determined up to N 
arbitrary phases entering B. The total number of inde- 
pendent variables used to parametrize in this way any 
density matrix p is equal to ./V 2 — 1. Since the literature 
seemed not to distinguish any natural measure in this 
space, we approached the problem by defining a product 
measure [Q 



Hu = Ai x u H . 



(2) 



The measure v is defined in the space of unitary matri- 
ces U(N), while A is defined in the (N — l)-dimensional 
simplex determined by the trace condition X)»=i di = 1. 
In Ref. [jjj we took for v the Haar measure on U(N), 
while the uniform measure Ai was used on the simplex. 
Our choice was motivated by the fact that the both com- 
ponent measures are rotationally invariant. For vr this 
follows directly from the definition of the Haar measure, 
while in Appendix A we prove that the uniform measure 
Ai corresponds to taking for the vector di the squared 
moduli of complex elements of a column or a row (say, 
the first column) of an auxiliary random unitary matrix 
V drawn with respect to vh 



In the sequel we will thus refer to the measure defined by 
(|2J) as the unitary product measure /z u . 

As correctly pointed out by Slater [@||, our choice of 
the measure is by far not the only one possible. He dis- 
cussed several possible measures, and proposed to pick 
the measure on the (N — 1)D simplex from a certain 
family of the Dirichlet distributions 
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where A > is a free parameter and C\ stands for a 
normalization constant. The last component is deter- 
mined by the trace condition d/v = 1 — d\ — . . . — djv— l- 
The uniform measure Ai corresponds to A = 1. Slater 
distinguishes also the case A = 1/2, which is related to 
the Fisher information metri c fl6|l , the Mahalonobis dis- 
tance |L7fl an d Jeffreys' prio r fig |, and was used for many 
years in different contexts |19J-21||. Since this measure 
is induced by squared elements of a column (a row) of 
a random orthogonal matrix (see Appendix A) , we shall 
refer to 
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as to the orthogonal product measure in the space of the 
mixed quantum states. Therefore both measures may 
be directly linked to the well known Gaussian unitary 
(orthogonal) ensembles of random matrices p2| , referred 
to as GUE and GOE. The measure p, u is determined by 
squared components of an eigenvector of a GUE matrix, 
while the measure /i may defined by components of an 
eigenvector of GOE matrices p3fl. Some properties of 
the orthogonal measure p have been recently studied in 
p4j . Let us stress that the name of the product measure 
(orthogonal or unitary) is related to the distribution A on 

the simplex d, while the random rotations U are always 
assumed to be distributed according to the Haar measure 
v H in U(N). 

It is interesting to consider the limiting cases of the 
distribution ([|). For A — v one obtains a singular dis- 
tribution concentrated on the pure states only km, while 
in the opposite limit A — > oo, the distribution is peaked 
on the maximally mixed state p* described by the vector 
d = {1/-/V, . . . , 1/-/V}. Changing the continuous param- 
eter A one can thus control the average purity of the 
generated mixed states. 



B. Separable states 

Consider a composite quantum system described by 
the density matrix p in the N dimensional Hilbert space 
H = Ha ® Hb . Dimension of the system N is equal to 
the product UatIb of the dimensions of the both subsys- 
tems. If the state p &H can be expressed as p — pa <8> pb , 



with pa € Ha and ps 6 TLb, it is called product 
state (or factorizable state). This occurs if and only if 
p = TibP^TtaP, where Tr^ and Tr# denote the oper- 
ations of partial tracing. In other words, for such states 
the description of the composite state is equivalent to the 
description in the both subsystems. 

A given quantum state p is called separable, if it can 
be represented by a sum of product states [J25| 



= Hp- 



QAi <X> QBi 



(6) 



where QAi and QBi are the states on TCb and Hb respec- 
tively. The smallest number k of the product states used 
in the above decomposition is called the cardinality of 
the separable state p [£6| . 

In general no explicit necessary and sufficient condi- 
tions are known for a mixed state to be separable. How- 
ever, Peres found a necessary condition showing that 
each separable state has the positive partial transpose 
0. Later Horodeccy demonstrated that for TV = 4 and 
TV = 6 this is also a sufficient condition |3| . To represent 
any state p it is convenient to use an arbitrary orthonor- 
mal product basis \ej) ® |e;), j = 1, . . . , ha; I = 1, • • ■ , Ub 
and to define the matrix Pjijw — (ej\ <S> (ez|p|e'-) (g> |e£). 
The operation of partial transposition is then defined B] 
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Even though the matrix p T2 depends on the particular 
based used, its eigenvalues {d[ > d' 2 >,... > d' N } do 
not. The matrix p T2 is positive iff all eigenvalues d! i are 
not negative. The practical application of the partial 
transpose criterion is thus straightforward: for a given 
state p one computes p T2 , diagonalizes it, and checks the 
signs of all eigenvalues. To characterize quantitatively 
the violation of positivity we introduced [lj] the negativity 
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which is equal to zero for all the states with positive 
partial transpose. 



C. Relative volume in the space of the density 
matrices 



the orthogonal and unitary product measures, and veri- 
fied, whether their partial transpose (0) are positive. The 
results are displayed in Fig.l as a function of the system 
size TV. Note that for TV > 6 we obtained in this way 
the volume Pt of the states with positive partial trans- 
position, which gives an upper bound for the volume of 
separable states. In fact Pt = Ps+Pb, where the volume 
Pb of the entangled states with positive partial transpose 
is studied in the section IV. 

The symbols are labeled by the size of the first subsys- 
tem riA- For both measures the symbols seem to lay 
on one curve, which would imply that Pt(^a,^b) — 
Pr(nA *tib)- However, this relation is only approximate, 
since Pt(2 X 6) ^ Pr(3 x 4) as pointed out by Smolin 
P7J . Numerical results for Pt and (t) for TV < 12 are col- 
lected in Table 1. The difference between Pt(2 x 6) and 
Pt(3 x 4) is not large, and was smaller than the statistical 
error of the results reported in W . Therefore it is reason- 
able to neglect for a while these subtle effects, depending 
on the way the TV dimensional system is composed, and 
to ask, how, in a first approximation, Pt changes with 
TV. 
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FIG. 1. Probability Pt of finding a state with positive par- 
tial transpose as a function of the dimension of the problem TV 
for unitary product measure (open symbols) and for orthog- 
onal product measure (full symbols). For TV < 6 it equals to 
the probability Ps of finding separable state, while for TV > 6 
it gives an upper bound for this quantity. Different symbols 
distinguish different sizes of one subsystem; ua = 2(o), 3(A) 
and 4(D). 



In M we presented several analytical lower and upper 
bounds for the volume of separable states. They were 
obtained assuming the unitary product measure, but the 
same reasoning can be repeated for other measures. The 
key result: an analytical proof that the volume of sep- 
arable states is positive and less than one, is obviously 
valid for any nonsingular measure. 

To analyze the influence of the measure chosen for the 
volume of the separable states P s we picked several ran- 
dom density matrices (circa 10 6 ) distributed according to 



Table 1. Probability Pt of finding a mixed states 
of size TV with positive partial transpose and the mean 
negativity (t) for two product measures orthogonal p, 
and unitary p u . For TV = 4 and TV = 6 one has Pt = Ps- 
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Figure 1, produced in a semilogarithmical scale, shows 
that for both measures the probability Pt decreases ex- 
ponentially with the system size N. Obtained numerical 
results allow us to conjecture that limjv^oo Pt(N) = 
for any (non-singular) probability measure used. Observe 
different slopes of both lines received for different proba- 
bility measures. The best fit gives Pt u ~ 1.8e _026Ar for 
the unitary product measure /x„, and Pto ~ 3.0e~ 055JV 
for the orthogonal product measure fi . Dependence of 
the probability Pt on the measure chosen is due to the 
fact that each measure distinguishes states of a different 
purity. This issue is discussed in details in the following 
sections. 



III. 2x2 CASE: POSITIVE PARTIAL 
TRANSPOSE ASSURES SEPARABILITY 

A. Purity versus separability 

For the TV = 4 case the partial transpose criterion 
is sufficient to assure the separability |||, so Pb = 
and Ps — Pt- Let us investigate, how the probabil- 
ity of drawing a separable state changes with its purity, 
which may be characterized by the von Neumann en- 
tropy Hi(g) = — Tr(ghxg). Another quantity, called the 
participation ratio 



R(g) = 



1 



Tr(p 2 )' 



(9) 



is often more convenient for calculations. It varies from 
unity (for pure states) to N (for the totally mixed state 
p* proportional to the identity matrix) and may be in- 
terpreted as an effective number of the states in the mix- 
ture. This quantity gives an lower bound for the rank 
r of the matrix p, namely r > R. Moreover, it is re- 
lated to the von Neumann-Renyi entropy of order two, 
Hiio) = InR(g). The latter, also called the purity of the 
state, together with other quantum Renyi entropies 
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is used, for q ^ 1, as a measure of how much a given state 
is mixed (see e.g. p8[). 

Figure 2 presents the probability distributions P(R) 
for N = 4 density matrices generated according to the 



both product measures. As discussed before, the or- 
thogonal measure [i is concentrated at less mixed states 
(lower values of R) then the unitary measure [i u . For 
example the mean value averaged over the orthogonal 
product measure (R) a ~ 2.184 is much smaller than the 
corresponding mean with respect to the unitary measure 
(R) u w 2.653. Observe a non-smooth behavior of both 
distributions at R = 3 (R = 2), for which the manifolds 
of a constant R start to touch the faces (edges) of the 3D 
simplex formed by d\ , di and d% . 
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FIG. 2. Purity and separability in iV=4-dimensional 
Hilbert space; open symbols represent averaging over the or- 
thogonal product measure \i while closed symbols are ob- 
tained with the unitary measure \i u \ a) probability distribu- 
tions P(R), b) conditional probability of finding a separable 
state as a function of the participation ratio R. All states 
beyond the dashed vertical line placed ai R — N ~ 1 = 3 are 
separable. 

Although the distributions P(R) differ considerably for 
both measures, the conditional probability of encounter- 
ing the separable state Ps(R) is almost measure indepen- 
dent, as shown in Fig. 2b. This is the main result of this 
section: the different results obtained for the probabil- 
ity Ps using various product measures fi\ are due to the 
different weights attributed to the mixed states. Since 



the average mixture (i?) MA grows monotonically with the 
parameter A (from 1 for A — > to 4 for A — > oo) also the 
probability P$ increases with this parameter from zero 
to unity. Note that for both curves the probability Ps 
achieves unity at R — 3: all sufficiently mixed states are 
separable. This fact has already been proved in Q, but 
see also [gj| for complementary, constructive results. 

Above considerations allow us to sketch the set of en- 
tangled states in the case TV = 4. In the analogy to the 
Bloch sphere, corresponding to TV — 2, we take the lib- 
erty to depict the set of all quantum states by a ball. 
Since it is hardly possible to draw a picture precisely 
representing the complex structure of the 15 dimensional 
space of the density matrices, Fig. 3 should be treated 
with a pinch of salt. In particular, the structure of the 
set of density matrices is not as simple and there exist 
several points inside the ball which do not correspond to 
density matrices. Furthermore, the 6-dimensional space 
of the pure states possesses the structure of the complex 
projective space CP 3 , much more complicated than a 
hypersphere. In the sense of the Hilbert - Schmidt met- 
ric, (A H s0i,/9 2 ) = y/T^[(pi - P2) 2 }), the set of the pure 
states forms a 6-dimensional subset of the 14-dimensional 
hypersphere of a radius \/3~/2 centered at p* = 1/4. 
Keeping this fact in mind, we represent this manifold 
by a circle in our oversimplified two dimensional sketch. 



B. Entanglement of formation 

After discussing the problem, how the probability of 
encountering a separable state changes with the degree 
of mixing R, we may discuss a related issue, how the av- 
erage entanglement depends on R. For this purpose we 
need a quantitative measure of entanglement of a given 
mixed state. Several such quantities have recently been 
proposed and analyzed |12||3(j-[40| , and none of them can 
be considered as the unique, canonical measure. How- 
ever, the quantity called entanglement of formation |l2[| 
plays an important role due to a simple interpretation: 
it gives a minimal amount of entanglement necessary to 
create a given density matrix. 

For a pure state \ip) one defines the von Neuman en- 
tropy of the reduced state 



E(ip) = -Trp A \np A = -Trp B lnp B , 



(11) 



where pa is the partial trace of \ip)(ip\ over the subsystem 
B, while ps has the analogous meaning. This quantity 
vanish for a product state. The entanglement of forma- 
tion of the mixed state p is then defined |12] 



E(p) =mmJ2 Pl E(^i) 



(12) 




pure states 
R=l 



and the minimum is taken over all possible decomposi- 
tions of the mixed state p into pure states 



FIG. 3. Sketch of the set of mixed quantum states for 
TV = 4: gray color represents the separable states 

The set of the separable states is visualized in Fig. 
3 as a 'needle of a compass': it is convex, has a positive 
measure and includes the vicinity of the maximally mixed 
state p» . Moreover, it touches the manifold of pure states 
(pure separable states do exist), but the measure of this 
common set is equal to zero. The more mixed state, 
(localized closer to the center of the 'ball'), the larger 
probability of encountering a separable state. All states 
with R > 3 are separable; for simplicity this complicated 
set is symbolically represented by a smaller circle. 
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(13) 



The decomposition of p into the smallest possible number 
of k pure states, for which this minimum is achieved, 
will be called optimal decomposition, while the number 
k will be called the cardinality of an entangled state. 
This definition may be considered as an extension of the 
concept of the cardinality of separable states introduced 
in |26|| , since for any separable state ps one has E(ps) = 
0. 

In Appendix B we present an algorithm allowing one to 
perform the minimization crucial for the definition (|12j). 
It gives an upper estimation of the entanglement of for- 
mation for an arbitrary density matrix of size TV. The 
algorithm proposed works fine for TV of the order of 10 or 
smaller. In the case of two qubits, discussed in this sec- 
tion, an analytical solution was found by Hill and Woof- 
ers fl||l5|] , who introduced the concept of concurrence. 

For any 4x4 density matrix p one defines the flipped 
state p — Op*0 T , where p* denotes the complex con- 
jugation and the orthogonal flipping matrix O contains 
only four nonzero elements along the antidiagonal: On = 
On = 1 and O23 = 032 = — 1- Concurrence C{p) is then 
defined ill 



C(p) :— max{0, ol\ — 0.2 — a^ — 0:4}, 



(14) 



where cms are the eigenvalues, in decreasing order, of the 
Hermitian matrix \TJppUp- Note that this matrix de- 
termines the Bures distance uM between p and p. In 
other words ens, are the non-negative square roots of the 
moduli of the complex eigenvalues of the non Hermitian 
matrix pp. 

Concurrence C of a given state p determines its entan- 
glement of formation |13| , |l5| 

E(p)=h(±[l + y/l-C*(p)]), (15) 

where 

h(x) :=-xln(x) - (l-x)ln(l-a;) (16) 

is the Shannon entropy of the 2-elements partition {x, 1 — 
x}. 
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FIG. 4. 2 x 2 system; a) The distributions P(E) of the 
entanglement of formation obtained for the density matrices 
generated according to p (white histogram, open symbols) 
and p u (gray histogram, closed symbols) , and the rotationally 
uniform distribution in the set of pure states (*); b) average 
entanglement E(R) (squares) and average negativity t(R) (di- 
amonds) for both measures p. 



Note that in the definition of entropy ( |ll| ) the natu- 
ral logarithm was used, (in contrast to the binary log- 
arithm present in the paper J13[|), so the entanglement 



E E [0,ln2]. Two histograms in Figure 4 present the 
probability distribution P(E) obtained for N = 4 ran- 
dom density matrices distributed according to both prod- 
uct measures p and p u . The singular peak at E — 0, cor- 
responding to the separable states, is omitted. Large en- 
tanglements of formation are rather unlikely. The mean 
values are not large: (E) sa 0.055 and (E) u « 0.018, 
since the averages are influenced by a considerable frac- 
tion of separable states with E = 0. The probability of 
obtaining a given value of E is larger for the orthogonal 
measure, which gives favor to more pure, and more likely 
entangled states. 

Both histograms may be compared with the probabil- 
ity distribution P(E) obtained for the ensemble of pure 
states, represented by stars in Fig. 4a. This distribu- 
tion is less peaked; vaguely speaking, different degrees of 
entanglement are almost equally likely among the pure 
states. The minimum of probability can be observed for 
maximally entangled states (E — In 2), while the mean 
(-E)purc ~ 0.328 is close to (ln2)/2. Since the singular 
distribution concentrated exclusively on pure states cor- 
responds to the case A — > in the distribution (0), we 
observe that the mean entanglement (E) decreases with 
increase of the parameter A, as the distributions A^ in- 
creasingly favor more mixed states. 

Although the mean entanglement (E) strongly de- 
pends on the measure used, the conditional mean en- 
tanglement E(R), averaged over all states of the same 
degree of mixing R, is not sensitive to the choice of a 
measure, as demonstrated in figure 4b. This allows us 
to formulate a general quantitative conclusion, valid for 
nonsingular measures in the space of density matrices: 
the larger the average degree of mixing R, the smaller 
the mean entanglement of formation E. For R > 3 one 
hasE(R) =0 |f. 



C. Negativity and concurrence 

In Ref. |jj we proposed a simple quantity t defined by 
(g), which characterizes quantitatively, to what extent 
the positivity of partial transpose is violated. As shown 
in Fig. 4b the conditional average t(R) does not depend 
on the measure applied and decreases monotonically with 
R. This dependence resembles the function E(R), which 
suggests a possible link between the both quantities. 

To analyze such a relation between these measures of 
entanglement, following the strategy of Eisert and Plenio 
p6[ , we generated 10 5 random density matrices p com- 
puting their concurrence C, entanglement E = E(C) and 
negativity t. As expected, the points at the plot E ver- 
sus t do not form a single curve. It means that both 
quantities, entanglement of formation and the negativ- 
ity, do not generate the same ordering in the space of 
4x4 density matrices. However, large correlation coef- 
ficients (approximately 0.978 for the orthogonal measure 
and 0.967 for the unitary measure) reveals a statistical 



connection between these measures. 
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FIG. 5. Ten thousands random density matrices of the size 
TV = 4 distributed according to the orthogonal product mea- 
sure; a) plot in the plane negativity - concurrence; b) plot of 
the difference C — t versus the participation R. 

It is particularly useful to look at the plane concurrence 
versus negativity. The data presented in Fig. 5a are ob- 
tained with the measure p . We observed, independently 
of the measure used, that all points are localized at or 
above the diagonal. This allows us to conjecture that for 
any density matrix p the following inequality holds 

t{p) < C{ P ). (17) 

A similar observation was already reported in pq| , where 



a modulus of the negative eigenvalue Em of the partially 
transposed matrix p T2 was used. Since for N = 4 no 
more than one eigenvalue d' 4 is negative |2q ], both quan- 
tities are equivalent and t — 2E^. Note that due to 
the conjecture ( |l7J ) we can attribute to negativity a more 
specific meaning. By means of equation ( |l5| ) and the fact 
that h(x) decreases for x > 1/2, negativity t allows us to 
obtain a lower bound for the entanglement of formation 
E. 

Numerical investigations show that the difference C — t 
is largest for mixed states with R ss 2, while it vanish for 
R > 3 and R = 1 (see Fig. 5b). In the former case 
all states are separable and C — t — 0. The latter case 
corresponds to pure states for which a-i = a^ = a^ = 



|5J and C = ai 



-2d' 



t. Thus the inequality (h 



becomes sharp for separable states or pure states. 



D. Mixed states with the same partition ratio R 

As demonstrated in Fig. 2b the conditional probabil- 
ity Ps (R) of encountering a separable state is similar for 
states with the same participation R, averaged over both 
product measures p and p u . This does not mean, how- 
ever, that the probability Ps is constant for each family 
of states p = UdU' defined by a given vector d with fixed 
participation ratio R. To illustrate this issue we discuss 
the case R = 2. 

Consider a vector of eigenvalues d with r nonzero ele- 
ments. This natural number (r e [1,4]) is just the rank 
of the matrix p. Any state p = UdU' can be expressed by 
the sum of r terms, pij = 5^j=i diUuU^. Moreover, the 
number of nonzero eigenvalues ai entering the definition 
of concurrence ( |14| ) equals to r |15|. 

Take any vector with r = 2 nonzero elements. In this 
case the formula (jlj) reduces to C = a± — a?. Since 
per definition a.\ > a.2, the concurrence is positive unless 
a\ = ai- Such degenerate cases occur with probabil- 
ity zero, (e.g. for diagonal rotation matrices U), so one 
arrives with a simple conclusion: 

For any set d of eigenvalues with r < 2 the probability 
Ps that a random state UdU^ is separable, is equal to 
zero. 

For concreteness consider three vectors of eigen- 
values characterized by r = 2,3 and 4. We put 
d a = {1/2,1/2,0,0}, d b = {2/3,1/6,1/6,0} and d c - 
{#1,3:2, #2, £2}) where X\ = (1 + \/3)/4 and xi — (1 — 
xi)/3. Each such vector generates an ensemble of density 
matrices p — UdU^ , where U stands for a random unitary 
rotation matrix of the size N — 4. Although all three en- 
sembles are characterized by the same participation ratio 
R — 2, the probabilities of generating an separable state 
are different. The case d a is characterized by r — 2, so 
Ps = 0. Numerical results obtained of a sample of 10 5 
random unitary matrices give Ps ~ 0.105 and 0.200 for 
df, and d c , respectively. Thus the probability Ps grows 
with the number r of pure states necessary to construct 



given mixed state p, or with the von Neuman entropy 

Hi. 

On the other hand, the average quantities character- 
izing entanglement (negativity, concurrence or entangle- 
ment of formation) decrease with r, provided, the partic- 
ipation R is fixed. For example, the mean entanglement, 
(E), equals 0.063, 0.057 and 0.042 for the ensembles d a , 
db and d c , respectively. Interestingly, in the latter case 
(or any other ensemble with g?2 = <^3 = di), one has 
«3 = a.A and C = t. 



IV. 2 x 4 CASE: POSITIVE PARTIAL 

TRANSPOSE DOES NOT ASSURE 

SEPARABILITY 

A. Purity and positive partial transpose 

For any system size the probability of finding the states 
with positive partial transpose depends on the measure 
used, as shown in Fig.l and Table 1. On the other hand, 
for any N the relations between purity and entanglement 
depend only weakly on the kind of the product measure 
used. 
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FIG. 6. As in Fig.2 for the 2 x 4 system (N = 8). The cir- 
cles in Fig. 6b represent the conditional probability of finding 
a state with positive partial transpose, Pt(R). Diamonds rep- 
resent the average negativity t(R) obtained with the measures 
p D (open symbols) and fi u (full symbols). 

Figure 6b presents the conditional probability of 
Pt{R) and the mean negativity t as a function of the 
degree if mixing R for the 2x4 system. For both quan- 
tities the results obtained with orthogonal and unitary 
product measures are difficult to distinguish. Thus the 
dependence of the total probability Pt on the measure 
used is strongly influenced by the likelihood of generating 
highly mixed states, described by the distribution P(R). 

These distributions for N — 8 are shown in Fig. 6a. 
The histogram for the unitary measure )jb u is shifted to 
larger values of R, with respect to the data obtained with 
the orthogonal measure. Quantitatively, the mean values 
read (R) u « 4.74 > (R) w 3.66. It is known that P T = 1 
for N > R — 1 |U . The right histogram, corresponding 
to the unitary measure, has a larger overlap with this 
region, which causes (Pt)u > (Pt)o- 

This observation is valid for an arbitrary matrix 
size, since for large N one has (R(N)) U rs N/2 while 
(R(N)) w N/3 {y|. For TV large enough the distribu- 
tions P(R) tend to Gaussians. They are centered at the 
mean values, which do depend on the measure, while the 
variance a 2 is of the order of N/5 for both measures under 
the consideration. Hence, the overlap with the interval 
[N — 1 , N] is larger for the measure p, u characterized by 
a larger mean value (R) u . 



B. Entanglement of formation 

Since for N > 4 there exist no analytical methods to 
compute the entanglement of formation of an arbitrary 
mixed state p, we relied on numerical computations. To 
perform the minimization present in the definition (fi~2) 
we worked out an algorithm based on a random walk in 
the space of unitary matrices U(M) with M > N. It is 
described in detail in Appendix B. Each run ends with 
an approximate optimal decomposition of the state p and 
provides an upper estimation of the entanglement E. To 
verify the accuracy of this technique we started with for 
the case N — 4, in which the explicit formula (fLa) is 
known. Computing numerically entanglement for 1000 
randomly chosen N = 4 mixed states we obtained the 
mean error of the order of 10~ 7 , while the maximal error 
was smaller than 10 -4 . 

At the beginning of each computation one has to 
choose the number M, determining the number of pure 
states in the decomposition. Since for N — 4 it is known 
that the cardinality of any state is not larger than 4 
Pq,E6[, it is sufficient to look for the optimal decompo- 
sition in the M = N = k = 4 dimensional space. For 
larger systems the problem of finding the maximal pos- 
sible cardinality is open. For each randomly generated 



mixed state p in the discussed 2x4 case we started to 
look for the optimal decomposition with At = N = 8, 
recorded the minimal entanglement Em=8 , an d repeated 
computations with M = 9, 10, . . . , M max . It is known 
[}42|,|7| that the maximal number of pure states does not 
exceed N 2 , but in practice we analyzed M £ [N, 2N] 

The number of degrees of freedom grows as M 2 , so the 
process of search for the optimal decomposition becomes 
less efficient with increase of the number M. However, for 
certain states we found better estimations for entangle- 
ment, e.g. EM=g(p) < Em=s{p)- I n these rare cases, the 
improvements of the estimations of E were very small, 
and repeating several times our procedure with M = 8 
the same upper bounds for entanglement of formation 
were reproduced. 

Thus our results do not contradict an appealing con- 
jecture that the cardinality k of any 2x4 mixed system 
is not larger than N = 8. Further work is still needed to 
verify, whether this conjecture is true. 

Note that the numerical algorithm to search for the op- 
timal decomposition and the entropy of formation, may 
also be used to look for the generalized entropy of for- 
mation E q , in the analogy to (|o|) and (|^), see |37J]. We 
found it interesting to study the quantity E2, which has 
a similar interpretation as the participation ratio R, and 
equals to unity for the separable states. 



estimations of the entanglement larger than E c are con- 
sidered as entangled) , the dependence of the results ob- 
tained on the cut-off value E c is weak. Moreover, these 
results do not depend on the exact values of the param- 
eters characterizing the random walk (see Appendix B). 
Consequently, we obtained an estimate of the volume of 
separable states for this case, Ps = Pt — Pb ~ 17.5%, as 
shown in the inset of Fig. 7. 



D. Bound entanglement and purity 

It is interesting to ask, whether a certain degree of mix- 
ing favors the probability of finding the bound entangled 
states. Grouping all 10 5 analyzed states in 30 bins ac- 
cording to the participation ratio R, we computed the 
conditional probabilities of entanglement. These results 
are shown in Fig. 7. Probability Ps increases mono- 
tonically with R while the probability of finding a free 
entangled state Pp = 1 — Pt decays with the partici- 
pation. On the other hand, the conditional probability 
Pb(R) of finding a bound entangled state exhibits a clear 
maximum at R ~ 5.5. If the mean purity is concerned, 
the bound entangled states are thus sandwiched between 
free entangled states (generally of high purity) and the 
separable states characterized by a high degree of mixing. 



C. Volume of the bound entangled states 

It is known H that for N — 8 there exist bound en- 
tangled states, which cannot be brought into the singlet 
form. All entangled states satisfying the partial trans- 
position criterion are bound entangled [J7|jqjl0||. It was 
shown in M that they occupy a positive volume Pg. 
Therefore Ps = Pt — Pb is smaller than the volume 
Pt of the states with positive partial transpose. Strictly 
speaking, the volume Pb of entangled states with positive 
partial transpose should be considered as a lower bound 
of the volume of bound entangled states, since it is not 
proven yet that all states with negative partial transpose 
are free entangled. 

To estimate Pb we generated 10 5 random density ma- 
trices of size N — 8. We worked with the unitary prod- 
uct measure /j, u , since, as shown in Table 1, the 2x4 
states chosen according to the orthogonal measure /i 
very seldom satisfy the partial transposition criterion. 
To save the computing time we estimated the entangle- 
ment of formation E only in the 2223 cases with posi- 
tive partial transpose. Setting an entanglement cut-off 
E c = 0.0003, (see appendix B), we found that 473 states 
enjoyed the entanglement E > E c . This gives a fraction 
of Pb « 4.7% of all states, or P B /P T = 21.3% of the 
states with positive partial transpose. Although these 
numbers are influenced by systematic errors (bound en- 
tangled states with E < E c are regarded as separable, 
while separable states with numerically obtained upper 
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FIG. 7. Conditional probabilities of finding the separable 
states (o), free entangled states (A) and bound entangled 
states (□) as a function of the participation ratio R. Results 
are obtained with 10 random density matrices of the size 
N — 8 distributed according to the measure fi u - The lines 
are drawn to guide the eye. The inset shows the pie-chart of 
the total probability of encountering separable states, bound 
entangled states (lower bound) , and free entangled states (up- 
per bound). 

Above results suggest that for the bound entangled 
states there exist a minimal participation ratio R or a 
minimal rank r. Preparing a sketch analogous to Fig. 3 
for A = 8, one should put the bound entangled states 
close to the center of the figure, but outside the symbolic 



'needle of the compass', which represents the separable 
states. 

It is worth noting, that the entanglement of formation 
for bound entangled states is rather small in compari- 
son to the mean entanglement of formation for the free 
entangled states, which violate the partial transposition 
criterion. The average taken over all free entangled states 
is (E)p ~ 0.05, the average taken over all bound entan- 
gled states is (E)b ~ 0.0033, (which is much larger than 
the cut-off value E c ). The maximal entanglement found 
for a bound state was only 0.0746. 



V. CLOSING REMARKS 

In this work we provide an attempt to characterize 
the statistical properties of the set of separable mixed 
quantum states. Certain level of caution is always rec- 
ommended by interpretation of any results of probabilis- 
tic calculations, especially if the space of the outcomes 
is infinite. Let us mention here the famous Bertrand's 
paradox: what is the probability that a randomly chosen 
chord of a circle is longer than the side of the equilateral 
triangle inscribed within the circle? The answer depends 
on the construction of the randomly chosen chord, which 
determines the measure in the infinite space of the pos- 
sible outcomes. 

Asking a question on the probability that a randomly 
chosen mixed state is separable, one should also expect 
that the answer will depend on the measure used. This 
is indeed the case, as demonstrated in this work for two 
products measures, and also shown by Slater |6J for a 
different measure related to the monotone metrics [ ff3[ . 
We reach, therefore, a simple conclusion, rather intuitive 
for an experimental physicist: the probability of finding 
an entangled state depends on the way, the states are 
prepared, which determines the measure in the space of 
mixed quantum states. 

On the other hand, in this paper we provide arguments 
supporting the conjecture that some statistical properties 
of entangled states are universal and to a large extend 
do not depend (or depend rather weakly) on the measure 
used. Let us mention only the exponential decay of the 
volume of the set of the separable states with the size 
N of the problem or the important relation between the 
purity of mixed quantum states and the probability to 
find a separable state. 

Studying the simplest case iV = 4we have shown that 
the distribution of entanglement of formation is close to 
uniform in [0,ln2] for all pure states. The more mixed 
states, the larger peak at small values of entanglement, 
the larger probability of finding a separable state. We 
have shown that the negativity t, a naive measure of en- 
tanglement, provides a lower bound for the entanglement 
of formation. 

Analyzing the more sophisticated problem N = 8 we 
developed an efficient numerical algorithm to estimate 



the entanglement of formation of any mixed state. In this 
way we could differentiate between separable states and 
the bound entangled states. About 79% of N = 8 states 
satisfying the positive transposition criterion are separa- 
ble. This result is obtained for random states generated 
according to the unitary product measure in the space 
of N = 8 density matrices, but we expect to get com- 
parable results for other, non-singular measures. Mean 
entanglement of formation for the bound entangled states 
is much smaller than for the free entangled states. Rel- 
ative probability of finding a bound entangled state for 
the 2x4 systems is largest for moderately mixed systems, 
characterized by the participation ratio close to R = 5.5. 
Even though this paper follows the previous work [[jj, 
the list of open problems in this field is still very long. 
Let us collect here some of them related to this work, 
mentioning also these, already discussed in the literature. 

A. N = 4, (2 x 2 systems) 

(i) Check, whether the dependence of the conditional 
probability on the participation ratio, Ps(R), obtained 
for two product measures (see Fig. 2b) holds also for the 
measures based on the monotone metrics or for the 
product Bures measure Q,Q. 

(ii) Prove the relation between the concurrence and the 
negativity: C > t. 

(iii) Find max(C — t) as a function of the participation 
ratio R, (see Fig. 4b). 

(iv) Check whether the following conjecture is true: 

If R(di) = R{d 2 ) and #i(di) > #i(cfe) then P s (fii) > 
Ps{p2)- The von Neuman entropy H\ and the partici- 
pation ratio R measure the degree of mixing of a given 
vector d, while Ps denotes the probability that a random 
state pi = UdiW is separable. 

(v) Find iso-probability surfaces in the simplex 
{^1,^2,^3} such that Ps(d) =const. 

B. N = 6, (2 x 3 or 3 x 2 systems) 

(vi) Find a lower bound for the entanglement of for- 
mation E, (in the analogy to negativity t, which gives a 
lower bound for C and E in the case N = 4) . 

(vii) Find an explicit formula for E in this case. 

C. N = 8, (2 x 4 or 4 x 2 systems) 

(viii) Find necessary and sufficient conditions for a 
bound entangled (or separable) state. 

(ix) Find the maximal entanglement of formation E of 
a bound entangled state. 

(x) Check whether the rank of bound entangled states 
is bounded from below. 

(xi) Check whether all states violating the partial 
transpose criterion are free entangled. 

(xii) Check whether the cardinality of any state is not 
larger than 8. 

D. General questions 

(xiii) Verify, whether the optimal decomposition of a 
given mixed state into a sum of pure states leading to 
the entanglement of formation, E = E\, also gives the 



10 



minimum of the generalized entanglement of formation 
E q . 

(xiv) For what Na x Nb composed systems the cardi- 
nality k of any mixed state in the N dimensional Hilbcrt 
space is less or equal to N — A^iVe? 

(xv) Check whether the entanglement of formation is 
additive. 

Not all of the above problems are of the same impor- 
tance. As most relevant we regard the questions (i),(viii) 
and the last two general problems. Preliminary results of 
Slater suggest pq] that the relation Ps(-R) for monotone 
metrics is similar to this obtained here for product met- 
ric, at least for N — 4. Concerning the question (viii): for 
separable states with N = 8, some necessary conditions, 
stronger than the positive partial transpose, are known 
p7| , ff8[ , but the sufficient conditions assuring the separa- 
bility are still most welcome. The problem of additivity 
of entanglement of formation is present in the literature 
(see e.g. ]15|]). Performing numerical estimations of the 
entanglement of formation E for several states of 2 x Nb 
systems we have not found any cases violating the state- 
ments (xiv) and (x v) [E9| . Recent results of Lewenstein, 
Cirac and Karnas [Eg suggest that the answer for the 
problem (xiv) is negative for the systems 3 x Nb with 
Nb > 3, but do not contradict this statement for the 
2 x Nb composed systems. Further effort is required to 
establish whether in this case the answer for the problem 
(xiv) is positive. 
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APPENDIX A: ROTATIONALLY INVARIANT 
PRODUCT MEASURES 

In this appendix we show that a vector of a N- 
dimensional random orthogonal (unitary) matrix gener- 
ates the Dirichlet measure (Q) with A = 1/2 (A = 1) 
in the (N — 1) D simplex. Although these results seem 
not to be new, we have not found them in the literature 
in this form and prove them here for convenience of the 
reader, starting with the simplest case N = 2. 

Lemma 1. Let O be a N x N random orthogonal ma- 
trix distributed according to the Haar measure on O(N). 



Then the vector d, 



Oil 



1, . . . , N is distributed 



according to the statistical measure on the (N — 1) di- 
mensional simplex (Dirichlet measure with A = l/2j. 



Proof. Due to the rotational invariance of the Haar 
measure on 0(A) the vector On is distributed uniformly 
on the N — 1 dimensional sphere S 1 ^^ 1 . Thus J2i=i di — 
1. 

For N — 2 the vector |Oii| is distributed uniformly 
along the quarter of the circle of radius 1. There- 
fore x = cos0, where e [0, 7r/2), and P(<fr) = 2/ir. 
Hence P(x) = P(cj))d(j)/dx = 2/(irVl-x 2 ). Another 
substitution £ = x 2 gives the required result: P(£) = 

P(x)dx/d£, = l/(7rV£(l - 0)- 

To discuss the general, N dimensional case, it is 
convenient to introduce the polar angles and to rep- 
resent any point belonging to the (N — 1) D sphere 



as Xn 



COSO N ^ 2 , P 



sintfjv-2, 



whcr 



e p- 



= 1 



Sj=i x i- Uniform distribution of the points on the 
sphere is described by the volume element dfl — 
sin ~ 2 0jv-2^at-2 • • • smOidOidcj). Changing the polar 
variables into Cartesian we obtain P(p) ~ 1/cos#at_2 = 
1/yl — p 2 ■ The last change of the variables £, := xf 
for i = 1, . . . , N allows us to receive P(£i, . . . , £/v-i) ~ 

[£i£2-.-£Ar-i(l-£i -6 6v-i)r 1/2 , which gives 

the statistical measure A 1 / 2 defined in Eq. (|j). □. 

Geometric interpretation of this result is particularly 
convincing for N = 3. Then the vector On covers uni- 
formly the sphere S 2 , while \On\ is distributed uniformly 
in the first octant. The points {^1,^2,^3} = {Cii^,^} 
lay at the plane z = 1 — x — y. Their projection into 
the x — y plane gives the statistical measure on the 2D 
simplex, i.e. the triangle {(0, 0), (1, 0), (0, 1)}. 

Lemma 2. Let U be a N x N random unitary ma- 
trix distributed according to the Haar measure on U(N). 
Then the vector di — \Un\ 2 , i = 1, ... ,7V is distributed 
according to the uniform measure on the (JV — 1 ) dimen- 
sional simplex (Dirichlet measure with A = 1). 

Proof. We will use the Hurwitz parametrization of 
U(N) [23] ], based on the angles tpki G [0, 7r/2] with 
0<k<l<N— 1. Their distribution can be deter- 
mined by the relation ipki =arcsin£ fc i \ , where £& are 
the auxiliary independent random numbers distributed 
uniformly in [0, 1] (see Ref. p3fl). 

In the simplest case N — 2 the vector d reads 
\Uii\ 2 = {cos 2 </? i,sin 2 (p 01 } = {£1, 1 - £1} and the vari- 
able d\ — £1 is distributed uniformly in the interval 
[0,1] (one dimensional simplex). For N — 3 one ob- 
tains d = {cos 2 ifi\2, sin </?i2 cos 2 <£>oi , sin ^12 sin t/'oi} = 

{1 — £ 2 1^2 (1 — £1)1^2 £i}j which is distributed uni- 
formly in the simplex {(0,0), (1,0), (0, 1)}. 

In the general A-dimensional case we 

get d = {COS 2 ifi N ^ 2 ,N~l, Sin 2 ip N ^ 2 .N-lCOS 2 ifiN^.N-i, 

sin 2 (p N -2,N-i sin 2 ^jv-3,jv-i cos 2 i/?jv-4,iV-i, 
sin 2 (p N -2,N-i ■ ■ ■ sin 2 ipi, N -i cos 2 ^ ,iV-i, 
sin 2 <£>at_2,jv-i ■ • • sin 2 </?i,jv-i sin 2 </5o,at-i}- Using uni- 
formly distributed random variables this vector may 

be written as {1 - ^^.^^(l - #<T°). 



SJV-1 ?AT-2 l 1 ?AT-3 )> 



A/(N-1) A/2,, 

SjV-1 '"'S2 l 1 
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&),&-: 



lff_! C2 £i}- This vector is uniformly dis- 

tributed in the N — 1 dimensional simplex, as explicitly 
shown in Appendix A of Q|. □. 

The above lemmae allow one to generate random 
points distributed in the simplex according to the both 
measures using vectors of random orthogonal (unitary) 
matrices. They may be constructed according to the al- 
gorithms presented in Ref. [E3[. Alternatively, one may 
take a random matrix of Gaussian orthogonal (unitary) 
ensemble, diagonalize it, and use one of its eigenvectors as 
in (f|). Random matrices pertaining to GOE (GUE) are 
obtained as symmetric (Hermitian) matrices with all ele- 
ments given by independent random Gaussian variables. 
Several ensembles interpolating between GOE and GUE 
are known p2j . Statistics of eigenvectors during such a 
transition were studied e.g. in |50| ; while the transitions 
between circular ensembles of unitary matrices were an- 
alyzed in p]]j . 



APPENDIX B: ENTANGLEMENT OF 
FORMATION - A NUMERICAL ALGORITHM 

1. Generating random density matrix 

In oder to generate an N X N random density ma- 
trix we write p = UdU^ and use the product measure 
fi = A\ x vh ■ The vector of eigenvalues d, taken accord- 
ing to the Dirichlet measure (Q), can be obtained from 
unitary random matrices as shown in Appendix A. Uni- 
tary random rotation matrix U distributed according to 
the Haar measure vh is generated with the algorithm 
presented in p3[ . Random state p, generated according 
to a given product measure, may be decomposed into a 
mixture of N pure states determined by its eigenvectors 



A ! 



Ei^x^i 



(Bl) 



j=i 



Note that the pure states |^/j) are not normalized to 
unity, but its norms are given by the eigenvalues di. 
Expansion coefficients of each of these states are given 
by the elements of the random rotation matrix; |\£j) = 

\/di{Uu, U2i, . . • , UNi}- 

There exist many other possible decompositions of the 
state p into a mixture of M pure states, with M > N. 
Let V be a random unitary matrix of size M distributed 
according to the Haar measure on U(M). Let V denotes 
a rectangular matrix constructed of the N first columns 
of V. Any such M x N matrix allows one to write a 
legitimate decomposition p' of the same state p 



M 



x>><& 



(B2) 



j=i 



\<Pi) 



N 

E 

m— 1 



Viml^m); i = 1, 



,M. 



(B3) 



Unitarity of the rotation matrix V assures the correct 
normalization Trp' = ^2 i= i(4>i\4>i) = 1. 

Assume that the composite N dimensional quantum 
system consists of two subsystems of the size Na and Nb , 
such that N = NaNb- It is then convenient to represent 
any vector \<pi) (of a non zero norm pi — {(j)i\<j)i)) , by 
a complex Na x Nb matrix A^\ which contains all N 
elements of this vector. To describe the reduction of the 
state \4>i) into the second subsystem we define a Nb x Nb 
Hermitian matrix 



B (i) .- [^Wjt^W. 



(B4) 



Diagonalizing it numerically we find its eigenvalues 
1 = 1, Nb ■ Rescaling them by the norm of the state 



5,« 



Pt we get b\ l > = by'/ pi satisfying Y,iJ\b\ 
compute the entropy of this partition 



di) 



N L 



£ B (W):=-^ftf»lntf 



1. We 



(B5) 



i=i 



giving the von Neuman entropy of the reduced state. En- 
tanglement of the state p' with respect to the rotated de- 
composition (B2) is equal to the average entropy of the 
pure states involved 



M 



E{p') = Y J P l E B {\^)) 1 



(B6) 



-,M 



where J2i=iPi = l - 

The entanglement of formation E of the state p is 
then defined as a minimal value Eb{p'), where the mini- 
murnis taken over the set of the decompositions p' given 
by (B2), (compare with the definition (12)). Rotation 
matrix V for which the minimum is achieved is called 
the optimal. Our task is to find the optimal matrix 
in the space of M dimensional unitary matrices where 
M = N,N + 1,...,N 2 . 

We found it interesting to consider also the generalized 
entanglement 



M 



where 



E q {p') = Y,Pi E S<t>i))> 



, N B 

E q (\4> i )):=-—hx(Y,W ) } 9 )- 
q i=i 



(B7) 



(B8) 



where 



The standard quantity Es{p) is obtained in the limit 

lim^i^p). 
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2. Search for the optimal rotation matrix 

The search for the optimal rotation matrix V has to 
be performed in the M 2 dimensional space of unitary 
matrices. Starting with M — N one has to consider 
the 16 dimensional space in the simplest case of N = 4. 
To obtain accurate minimization results in such a large 
space one should try to perform some more sophisticated 
minimization schemes, for example the stimulated an- 
nealing. Fortunately, the optimal rotation matrix V is 
determined up to a diagonal unitary matrix containing 
M arbitrary phases. Therefore, one can hope to get rea- 
sonable results with a simple random walk, moving only 
then if the entanglement decreases. Performing only the 
'down' movements in the M 2 dimensional space one has a 
good chance to land close to the M dimensional manifold 
defined by optimal matrices equivalent to V . This corre- 
sponds to fixing the temperature to zero in the annealing 
scheme, and simplifies the search algorithm. 

To perform small movements in the space of unitary 
matrices we will use M x M Hermitian random matrices 
H, pertaining to the Gaussian Unitary Ensemble (GUE). 
They can be constructed by independent Gaussian vari- 
ables with zero mean and the variance (<r m c n ) 2 = (1 + 
S mn )/M for the real part and (cr^) 2 = (1 — S mn )/M for 
the imaginary part of each complex element H mn — H* m . 
We generate random matrix H and take W = e txIi as a 
unitary matrix, which might be arbitrarily close to the 
identity matrix. Our strategy consists in fixing the initial 
angle xo, performing random movements of this size and 
then gradually decreasing the angle %• 

The detailed algorithm of estimating the entanglement 
of formation of a given N x N state p is listed below. 

(1) Fix the nu mber M of the components of the de- 
composition (B2). Start with M = N. 

(2) Generate random unitary rotation matr ix V of size 
M, which defines the decomposition p' in ( |B2[). C ompute 
the entanglement E — Eb(p') according to (B5,B6). 

(3) Set the initial angle x — Xo- 

(4) Generate a random M X M GUE matrix H and 
compute V' — Uexp(ix-ff). Calculate the entanglement 
E' for the decomposition p' generated by V. 

(5) If E' < E accept the move (substitute V := V and 
E := E') and continue with the step (4). In the other 
case repeat the steps (4-5) I c hange times. 

(6) Decrease the angle x '■— a X, where a < 1. 

(7) Repeat the steps (4-6) until x < Xend- Memorize 
the final value od the entanglement E. 

(8) Repeat L mat times the steps (2-7) starting from a 
different initial random matrix V. 

(9) Memorize the value Em, defined as the smallest of 
Lmat repetitions of the above procedure. 

(10) Set M =: M + 1 and repeat the steps (2-9) until 
M = M max . 

(11) Find the smallest value of Em, M = N, . . . , M max . 
This value E m i n = Em, gives the upper bound for the 
entanglement of formation of the mixed state p, while the 



size M* of the optimal rotation V* may be considered as 
the cardinality of p. 



3. Remarks on estimating the entanglement of 
formation 

The accuracy of the above algorithm may be easily 
tested for the case N = 4, for which the analytical for- 
mula (pit) exist. Results mentioned in the section IV B, 
giving the mean error of the estimation of the entangle- 
ment smaller than 10 -7 , were obtained with the following 
algorithm parameters: the initial angle xa — 0.3, the fi- 
nal angle Xend = 0.0001, the angle reduction coefficient 
a = 2/3, the number of iterations with the angle fixed 
Ichange = 25, the number of realizations L mat = 3. Using 
relatively slow routines interpreted by Matlab on a stan- 
dard laptop computer we needed a couple of minutes to 
get the entanglement of any mixed state p. Although we 
performed test searches with M = 4, 5, . . . , 8 the optimal 
rotation was always found for M = N = 4. 

The same algorithm was used for random states with 
N = 2 x 4 = 8. In this case the simplest search with 
M = N = 8, performed in 64 dimensional space, requires 
much more computing time. It depends on all parame- 
ters characterizing the algorithm; one may impose there- 
fore an additional bound on the total number / of gen- 
erated random matrices V . To estimate the volume of 
the bound entangled states we performed the above algo- 
rithm only for the states with positive partial transpose. 
Setting the final angle at Xend = 0.0002 we obtained in 
a histogram P(E) a flat local minimum at E m ~ 0.0003. 
The minimum is located just right to the singular peak 
at even smaller values of E, corresponding to separable 
states. The cumulative distribution P c = J E P(E)dE 
was found not to be very sensitive on the position of the 
minimum E m . We could, therefore, set the cut-off value 
E c to the position of minimum E m , and interpret the 
quantity P c as the relative volume of the bound entangled 
states. In the computations described in section IV C we 



took xo = 0.3, a = 2/3, I cha 



25 and L, t 



5, and 



for M = N = 8 obtained the mean number of iterations 
(I) of the order 5 x 10 4 . 

The other possibility to distinguish the bound entan- 
gled states from the separable states consists in studying 
the dependence of the obtained upper bound on the en- 
tanglement E on the total number I of iterations per- 
formed. Numerical results obtained for the separable 
states show that E decreases with the computation time 
not slower than E(I) = a/I. Assuming a similar ef- 
fectiveness of the algorithm for the nonseparable states 
(with a non zero entanglement of formation Ef orm ), we 
have E(I) — Ej orm + b/I. This allows us to design a sim- 
ple auxiliary criterion: the state p is separable if (for all 
realizations of the random walk starting from the differ- 
ent matrices V) for sufficiently large number of iterations 
/ one has E(I) < E(I /2)/2. If this condition is not ful- 
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filled the state p can be regarded as entangled. Using 
this method we obtained the estimation for the volume 
of bound entangled states similar to P c . 
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